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Motivated by the recent experimental success in realizing synthetic spin-orbit coupling in ultracold 
atomic systems, we consider JV-component atoms coupled to a non-Abelian SU(iV) gauge field. 
More specifically, we focus on the case, referred to here as "SU(3) spin-orbit-coupling," where the 
internal states of three-component atoms are coupled to their momenta via a matrix structure that 
involves the Gell-Mann matrices (in contrast to the Pauli matrices in conventional SU(2) spin- 
orbit-coupled systems). It is shown that the SU(3) spin-orbit-coupling gives rise to qualitatively 
different phenomena and in particular we find that even a homogeneous SU(3) field on a simple 
square lattice enables a topologically non-trivial state to exist, while such SU(2) systems always 
have trivial topology. In deriving this result, we first establish an exact equivalence between the 
Hofstadter model with a 1/N Abelian flux per plaquette and a homogeneous SU(iV) non-Abelian 
model. The former is known to have a topological spectrum for N > 2, which is thus inherited by 
the latter. It is explicitly verified by an exact calculation for N = 3, where we develop and use a 
new algebraic method to calculate topological indices in the SU(3) case. Finally, we consider a strip 
geometry and establish the existence of three gapless edge states - the hallmark feature of such an 
SU(3) topological insulator. 



Introduction. Spin-orbit coupling has a long history 
in solid-state systems and is responsible for a variety 
of important topics of current interest. Notable exam- 
ples include applications in spintronics [1] and the real- 
ization of topological insulators [H [3]. Recently, inter- 
est in spin-orbit coupling has come to the fore in the 
seemingly disparate area of ultracold atoms with the ad- 
vent of synthetic gauge fields [3]. Such gauge fields have 
been employed to mimic magnetic fields [5-8 a as well as 
Rashba-Dresselhaus spin-orbit coupling in both bosons 
[5] and fermions [TU] . This progress opens doors not only 
to quantum simulation of spin-orbit coupled solid-state 
systems [TTJ [T^], but also to the realization of a much 
larger class of structures that can be engineered in the 
ultracold laboratory but do not exist in the solid state 
(see, for instance, [T3HT5] ). 

The conventional spin-orbit coupling manifests itself 
as a Zeeman magnetic field that depends on the elec- 
tron's momentum. Hence a typical spin-orbit term in a 
continuum model of a solid is b(p) • <r, where the form 
of the momentum-dependent internal field, b(p), is dic- 
tated by symmetries of the crystal structure and <x is a 
vector of Pauli matrices, which mathematically are gen- 
erators of the SU(2) group that act on the electron's 
SU(2) spin. In contrast to solids, synthetic spin-orbit 
structures in ultracold atoms are built from the ground 
up and are not constrained by fundamental symmetries. 
Furthermore, since the "spin" itself is synthetic, there is 
no requirement that it be associated with a representa- 
tion of the SU(2) group. Hence, a much larger space of 
SU(iV) spin-orbit couplings become available for multi- 
component atoms, ^2i = ib l (p)Xi, where Xi are in prin- 
ciple any of the (iV 2 — 1) Hermitian generators of SU(A r ) 
(e.g., the Pauli matrices for N = 2, the Gell-Mann ma- 
trices for N = 3, etc.). 



In this work we introduce the notion of SU(7V) spin- 
orbit coupling and also develop an algebraic formalism 
to explore topological structures in such lattice systems. 
Because the class of realizable SU(iV) gauge fields is ex- 
tremely large, we focus on a specific SU(iV > 2) model 
that has several attractive features: (i) It is qualitatively 
different from its SU(2) counterpart; (ii) It has a non- 
trivial topological structure with a non-zero Chern num- 
ber |16j ; and (iii) It is very simple and is potentially real- 
izable in experiment (note that thus far only spatially ho- 
mogeneous gauge fields have been realized experimentally 
[HIHD])- More specifically we consider a simple square lat- 
tice in two-dimensions with nearest-neighbor hoppings 
only and a homogeneous (i.e., position- independent) 
SU(iV) gauge field. The corresponding Bloch Hamilto- 



H(k) = —2t cos(k x — A x ) + cos(k y — A y 



(1) 



where t is the hopping and the gauge fields A Xjy are con- 
stant N x N Hermitian matrices. We demonstrate that 
there is a qualitative difference between such SU(2) and 
SU(3) gauge fields, where the latter have a much richer 
algebraic and geometric structure. In particular, for the 
familiar SU(2) case, it can be seen that there are no gauge 
fields which result in non-zero Chern numbers |17j . On 
the other hand, we show by explicit construction that 
SU(3) gauge fields with non-zero Chern numbers do ex- 
ist. 

To construct such an "SU(3) topological insulator" we 
relate the model ([I]) to the Hofstadter model [T5], which 
describes particles hopping on a square lattice under a 
uniform magnetic field (but non-uniform gauge field). 
This is achieved by removing the spatial dependence of 
the Abelian gauge field of the Hofstadter model by a 
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gauge "twisting" transformation. We then move on to 
compute the Berry curvature and Chern numbers for gen- 
eral SU(3) systems. In particular, we derive an equation 
for the Berry curvature which is directly expressed in 
terms of quantities of the Bloch Hamiltonian. This for- 
mula provides a generalization of a useful and well-known 
expression for SU(2) systems [c.f. Eq. {[sj ] . We apply this 
method to find the curvature and Chern numbers of the 
SU(3) model. To elucidate the behavior of the edge states 
of this model, we investigate the system in a strip geom- 
etry. Finally we comment on the experimental feasibility 
of SU(3) systems. 

Mapping from the Hofstadter Model. We will now give 
the details of the mapping of the Hofstadter Model to a 
Bloch Hamiltonian having the form of Eq. ([lj . The Hof- 
stadter Model [18] describes particles on a square lattice 
under an Abelian, but spatially-dependent gauge field 
describing a uniform magnetic field. To extend the Hofs- 
tadter model to SU(iV) systems we consider its N decou- 
pled replicas, each having the same flux per plaquette. 
Our starting point is thus the Hamiltonian 

Hum = -t + *J e -«™(**+^. + . + H . c .) , 

i 

(2) 

In this equation = (ipn, ipa., • ■ • , ipiN) T are SU(A) 
spinor operators, S z — diag (s,s — 1,...,— s) where 
2s + 1 = N, x and y are the two square lattice vec- 
tors where the lattice constant is set to unity, sc< = x • r.; 
where r< is the position of the i th lattice site, and a gives 
the magnitude of the flux. Since, S z is diagonal in this 
representation, the model trivially decouples into N in- 
dependent copies of the Hofstadter model, each having 
2na flux per plaquette. We restrict the flux to be related 
to the number of spin components as a — 1/N. 

We will first illustrate the mapping for the case of 
two-component spins and later describe how to gener- 
alize. For this case we apply the gauge transformation 
"fi — > e~ l % (7:cXil i>i where a x is a Pauli matrix. This 
transformation rotates the spinors about the x-axis by 
a position-dependent angle. As can be seen after some 
straightforward algebra, this transformation removes the 
spatial dependence of the second term in Eq. ([2| at the 
cost of introducing a non- Abelian x-component into the 
gauge field. In particular, after the gauge transformation 
Eq. |2| becomes 

H = -tJ2 (^le-^i+x + ^\e- lA ^. l+9 + H.c.) 

i 

= £*!#(k)*k, (3) 

k 

where in the second line we have taken the Fourier trans- 
form. For this case, the non-Abelian gauge fields of 
Eq. |lj can be expressed as Pauli matrices as (A x , A y ) = 
\(p x ,o z ). The Bloch Hamiltonian can also be expanded 



and rewritten as H(k) — — 2t[sin(k x )a x + sm(k y )a z ] 
which is a lattice version of Rashba spin-orbit coupling. 

We now generalize this mapping to any integer N. As 
before, we perform a gauge transformation ^ — > W 1 '*;. 
The unitary matrix U is defined to have = U n +i^ n = 
—i for 1 < n < N — 1 with zeroes elsewhere. One can 
verify for this matrix that e - i2 ™ S *U = e -«*-a(&,-i). 
Therefore this gauge transformation will completely re- 
move the position dependence of the second term in 
Eq. ([2]), and the transformed Hamiltonian will corre- 
spond to particles on a square lattice under a homoge- 
neous gauge field. 

The SU(3) case of this general mapping will be con- 
sidered in detail below. For this case, the non-Abelian 
gauge fields arrived at through the mapping which enters 
Eq. can be expressed in terms of Gell-Mann matrices 
PI as 

A x = ^(^2 - A 5 + %i) and A y = |(A 3 + V3X S ). (4) 

By expanding the cosines, the Bloch Hamilto- 
nian can also be written as 'H(k) = b(k) • A 
where b(k) = — t (cos(k x ), sin(fc z ), ^ sm(k y ) — 
\ cos(ky), cos(k x ), - sin(fc x ), cos(/c x ), sin(fc x ), ^ cos(fc !/ ) + 
I sin(fc y ) ) is an eight-component vector and A is a vector 
composed of the eight Gell-Mann 3x3 matrices. 

Geometrical Method for Berry Curvature and Chern 
Number Computation. We now describe a geometrical 
method of computing the Berry curvature and Chern 
numbers for general SU(3) systems. We first write down 
expressions which are valid for any N. The Berry cur- 
vature 0„(k) [20] is defined in terms of the normalized 
eigenstates Xkn of the Bloch Hamiltonian as 

n»(k) = i (d kx xt n dk y Xkn ~ <9 fcy Xk,A*Xk«) > (5) 

where n labels the eigenstate (or band). The Chern num- 
ber for a particular band is defined as [16] 

v n = ^ [ d 2 kn n (k), (6) 

*n Jbz 

where the integral is performed over the first Brillouin 
zone (BZ). The Berry curvature can also be expressed in 
terms of eigenstate projection operators P^ n — Xk« <8> Xk„ 
where ® denotes the outer product through the useful 
relation [2"Tl |2"2] 

fl n (k)dk x A dk y = iTr(P kn A dP kn A dP kn ), (7) 

where dk x A dk y — —dk y A dk x . 

Before generalizing we first describe a well-known ge- 
ometrical expression for the Berry curvature for SU(2) 
systems (see, e.g., [2]). This will be used to demonstrate 
that SU(2) Bloch Hamiltonians of the form Eq. are in 
general topologically trivial. The Bloch Hamiltonian for 
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SU(2) systems can be expressed in terms of Pauli matri- 
ces as 'H(k) = a(k) + b(k) • a. The projection operators 
corresponding to the two eigenstates can be written in 
terms of b(k) as P k ± = \ [1 ± b(k) • <x/|b(k)|] . Insert- 
ing this into Eq. ([7| then gives 



fi±(k) 



2|b(k)|- 



rb(k) • [9 fex b(k) X %b(k)] . (8) 



Thus, the Berry curvature can be expressed directly in 
terms of the Bloch Hamiltonian, rendering the interme- 
diate steps of computing its eigenstates and evaluating 
Eq. |5| unnecessary. For SU(2) systems, one can write 
arbitrary gauge fields of H(k) as linear combinations of 
Pauli matrices as A XjV — u x _ y + ~v x ,y ■ o - - After expanding 
the exponents to obtain b(k), it is a straightforward ex- 
ercise to verify that 9fc x b(k) x dk b(k) a v x x v y . Then 
through Eq. Q one sees that the Berry curvature van- 
ishes identically rendering SU(2) systems described by 
Eq. |l]) topologically trivial. 

We now move on to develop a central technical re- 
sult of our work, namely the generalization of Eq. |8]) 
to SU(3) systems. We will utilize the elegant formal- 
ism presented in |23| which describes an efficient way to 
represent pure-state density matrices (or projection op- 
erators) for three-state systems. For SU(3) systems, a 
general Bloch Hamiltonian can be expressed in terms of 
the eight Gell-Mann matrices as 



H(k) = o(k) + b(k) • A, 



(9) 



where a(k) is a scalar and b(k) is an eight-dimensional 
real vector. The product of two Gell-Mann matrices can 
be written as A Q A fc = |5 ab + d abc X c + if abc \ c where d abc 
and f a bc are the symmetric and antisymmetric structure 



constants of SU(3) [H]- These structure constants define 
three bilinear operations for the eight-component vec- 
tors. In particular, one has the dot product u- v = u a v a , 
the cross product (u x v) a = f a bc u bV c , and the so-called 
star product [33] (u * v) a = \f3d a bcUbV c for two arbi- 
trary vectors u and v where repeated indices are summed 
over. One can also write eigenstate projection operators 
in terms of the Gell-Mann matrices as 



Pk» = Xk« ® xt n = ^(1 + V3n kn • A) 



(10) 



where TrPk„ = 1. The condition that (Pknj = Pkn> 
leads to two constraints on the vector n^ n which are nk« • 
nkn = 1 and n^ n * rikn = nk n [23]. Due to the star- 
product constraint, nk« lies in a restricted region of S 7 . 
This can be compared to the SU(2) system where the 
vector analogous to n k „ can lie anywhere in S 2 . 

Now we will express n^ n in terms of b(k) appearing in 
the Bloch Hamiltonian Eq. For projection operators 
corresponding to eigenstates we have [J\„,H(k)] = so 
that b(k) x nkn = 0. One can verify that this equation, 
along with the above constraints, is satisfied by nk n = 
£kn [7knb(k) + b(k) * b(k)] with coefficients 

7k „ = 2|b(k)|cos^ k + y^ ; (11) 



£k« 



|b(k)|2[4cos 2 (0 k + f n)-l] 



where #k = | arccos 



b(k)-b(k)*b(k) 
|b(k)|3 



and n runs from one 



to three. The resulting expression for Pk„ can be inserted 
into Eq. Q to obtain the Berry curvature. One finds 



fi *( k ) = ~^J2 b 2d k* h x ^b + jd k b x d ky (b * b) + 7 9 fex (b * b) x d k b + d kx (b * b) x d ky (b * b)] ■ ( 7 b + b * b), 

(12) 



where we have suppressed the k, n arguments on the 
right-hand side. Notice that due to orthogonality rela- 
tions, the derivatives do not act on the coefficients. While 



Eq. (12 1 is complicated in appearance, it is straightfor- 
ward to compute with a given b(k). This equation pro- 
vides an explicit expression for the Berry curvature in 
terms of quantities from the Bloch Hamiltonian and thus 
should be viewed as a generalization of Eq. (JsJ) to SU(3) 
systems. 

Analysis of SU(3) model. Having established the above 
formalism, we now move on to analyze the specific SU(3) 
model arrived at above, given by Eqns. and Q. The 
resulting b(k) can be directly inserted into Eqns. (11) 



and ( 12 1 to find the Berry curvature for this system. One 
finds 



n„(k) 



2cos(46» k + ^fn) -3 
V3[l + 2cos(26» k - ^n)}' 



where #k = \ arccos 



(13) 



In addi- 



^| (cos^fca) + cos(3A^)) 

tion, using the expression E^ n — Tr{Pk n H(k)}, the bulk 
eigenenergies are found to be 



= k cos(3g k ) + 2cos(fl k + f n) 
k " l + 2cos(26» k + \n) ' 



(14) 
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FIG. 1. The spectrum of SU(3) model in a strip geometry. 
The bulk states correspond to black lines while the edge states 
correspond to the thick red lines. 



These bands are gapped and ordered such that £^1 < 
E\z2 < £-k3- With the above expressions for the curva- 
ture, the Chern numbers can be computed via Eq. Q 
and are found to be v^) = (—3, 6, —3). 

Due to the bulk-boundary correspondence [23], non- 
zero Chern numbers imply the existence of edge states. 
To elucidate the behavior of these edge states, we inves- 
tigate the SU(3) system in a strip geometry. We apply 
periodic boundary conditions in the y-direction, and take 
a system of finite length in the x-direction. The system 
in this strip geometry is described by 

ffstrip = ~t^2 [^H k v) 2cos ( k v - (15) 



— iA 



%+i(k y ) + *l +1 (k y )e iA *%(k y ) 



where i now is a one-dimensional finite sum. The eigen- 
states of -Hstrip are plotted in Fig. [T] The spectrum ex- 
hibits characteristic topological edge states that connect 
the bands with different Chern numbers. 

In conclusion we make a few general remarks. First, 
we note that while the SU(3) topological insulator con- 
structed here relies on spin-orbit coupling of a new type 
and while the calculation of Chern numbers requires a 
new algebraic construction, its overall topological char- 
acterization resides within the existing general classifica- 
tion scheme [25l [26] and corresponds there to a lattice 
quantum Hall state labelled by an integer topological in- 
dex. However in contrast to solid-state systems where 
the absence or presence of time-reversal symmetry is an 
obvious physical constraint, for synthetic spin-orbit sys- 
tems the notion of time-reversal symmetry does not have 
such a direct meaning, because the synthetic spins do not 
behave like real spins under time reversal. Classification 
of cold-atom Hamiltonians with respect to transforma- 



tions of time-reversal type can still be formulated but in 
a more formal way by examining the existence of an anti- 
unitary symmetry of the Hamiltonian which may or may 
not have a direct physical interpretation. From this per- 
spective, our Hamiltonian does not have such a symme- 
try. One can argue that in general such Chern topologi- 
cal insulators are much easier to realize with cold atoms 
than Zi topological insulators, because imposing an ad- 
ditional unphysical symmetry would require fine-tuning 
the synthetic Hamiltonian, in contrast to the situation in 
the solid state where in the absence of external magnetic 
fields and magnetic impurities time-reversal invariance 
is automatically preserved. Finally, we briefly comment 
on experimental realization and detection of the SU(3) 
system. There exists a considerable literature on the re- 
alization of synthetic gauge fields in cold atom systems 
(for a review, see [IS])- The gauge fields from Eq. Q 
can be realized with variations of the so-called iV-pod 
schemes. The experimental observation of the proposed 
state can be achieved by the same methods as were pro- 
posed to observe more conventional quantized Hall states 
(see, e.g., [13 |2g]). 

Acknowledgements - This research was supported by 
JQI-PFC (RB), ARO-MURI (GB), and US-ARO (VG). 



[1 

[2 
[3 

K 

[5] 
[6] 

[8 

[9 
[10 

[11 

[12 
[13 



I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011) . 

Y.-J. Lin, R. L. Compton, A. R. Perry, W. D. Phillips, 
J. V. Porto, and I. B. Spielman, Phys. Rev. Lett. 102, 
130401 (2009). 

Y. J. Lin, R. L. Compton, K. Jimenez-Garcfa, J. V. 
Porto, and I. B. Spielman, Nature 462, 628 (2009). 
M. Aidelsburger, M. Atala, S. Nascimbene, S. Trotzky, 
Y. A. Chen, and I. Bloch, Phys. Rev. Lett. 107, 255301 

(2011) . 

J. Struck, C. Olschlger, M. Weinberg, P. Hauke, J. Si- 
monet, A. Eckardtand, M. Lewenstein, K. Sengstock, 
and P. Windpassinger, Phys. Rev. Lett. 108, 225304 

(2012) . 

K. Jimenez-Garcfa, L. J. LeBlanc, R. A. Williams, M. C. 
Beeler, A. R. Perry, and I. B. Spielman, Phys. Rev. Lett. 
108, 225303 (2012). 

Y.-J. Lin, K. Jimenez-Garcfa, and I. B. Spielman, Nature 
471, 83 (2011). 

L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T. Yef- 
sah, W. S. Bakr, and M. W. Zwierlein, arXiv: 1205.3483 
(2012). 

L. Mazza, A. Bermudez, N. Goldman, M. Rizzi, M.-A. 
Martin-Delgado, and M. Lewenstein, New J. of Phys. 
14, 015007 (2012). 

I. Bloch, J. Dalibard, and S. Nascimbene, Nature Phys. 
8, 267 (2012). 

J. Dalibard, F. Gerbier, G. Juzeliunas, and P. Ohberg, 



5 



Rev. Mod. Phys. 83, 1523 (2011). 
[14] N. R. Cooper, Phys. Rev. Lett. 106, 175301 (2011). 
[15] B. M. Anderson, G. Juzeliunas, V. M. Galitski, , and 

I. B. Spielman, Phys. Rev. Lett. 108, 235301 (2012). 
[16] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 

M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 
[17] One should note, however, that by relaxing the condition 

of spatially-uniform gauge field, it is possible to obtain 

topological states [29] . 
[18] D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 
[19] H. Georgi, Lie Algebras In Particle Physics: from Isospin 

To Unified Theories (Westview Press, 1999). 
[20] M. V. Berry, Proc. R. Soc. London, Ser. A 392, 45 (1984). 
[21] J. E. Avron, R. Seiler, and B. Simon, Phys. Rev. Lett. 



51, 51 (1983). 

[22] J. E. Moore and L. Balents, Phys. Rev. B 75, 121306 
(2007). 

[23] G. Khanna, S. Mukhopadhyay, R. Simon, and 

N. Mukunda, Ann. Phys. 253, 55 (1997). 
[24] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 
[25] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. 

Ludwig, Phys. Rev. B 78, 195125 (2008). 
[26] A. Kitaev, AIP Conf. Proc. 1134 (2009). 
[27] T. D. Stanescu, V. Galitski, J. Y. Vaishnav, C. W. Clark, 

and S. Das Sarma, Phys. Rev. A 79, 053639 (2009). 
[28] E. Zhao, N. Bray-Ali, C. J. Williams, I. B. Spielman, 

and I. I. Satija, Phys. Rev. A 84, 063629 (2011). 
[29] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 



